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A.  ABSTRACT 
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Algorithms  were  developed  for  the  generation  of  grids  in 
2D  and  3D.  The  3D  code  development  involved  extensive  use  of 
computer  Symbolic  Manipulation.  A  rigorous  code  validation 
procedure  was  developed.  A  controversy  on  the  effect  of  strong 
coordinate  stretching  was  resolved,  with  proof  that  the  order  of 
accuracy  is  not  reduced.  The  contention  that  a  well-known 
method  always  produces  a  nan-folded  grid  was  disproven.  A 
family  of  adaptive  algorithms  was  developed,  involving  both 
interior  adaption  and  especially  boundary  adaption.  Initial 
work  was  begun  on  an  extension  of  the  variational  approach  to 
grid  generation.  The  grid  generation  algorithms  have  been 
applied  in  the  2D  and  3D  ELF  codes  (ELectric  Field)  used  in  the 
design  of  laser  electrodes  and  switches.  Applications  are  in 
the  Pulsed  Power  area,  including  SDI  research  and  development. 


B.  STATEMENT  OF  WORK 


"Compatible  with  the  level  of  funding,  the  contractor 
will  exert  his  best  efforts  to  develop  numerical  methods  for 
adaptive  grid  generation  for  two-  and  three-dimensional  internal 
and  external  flows.  The  numerical  methods  developed  under  the 
proposed  contract  will  be  documented  in  open— 1 i ter ature  publica¬ 


tions 


presented 


scientific  and  engineering  meetings. 


sample  codes  will  be  made  available 
government  agencies  and  contractors. " 


interested  users  within 


1 
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C.  STATUS  OF  THE  RESEARCH 

The  algorithms  and  code  for  the  generation  of  the  2D  grid 
generation  equations,  the  solution  of  the  finite  difference 
equations,  and  the  verification  process  (see  below)  were  all 
successfully  developed  using  hand  coding.  The  solution  methods 
included  hopscotch  SOR  and  the  much  more  efficient  (but  geometry 
sensitive)  marching/semidirect  methods. 


The  3D  algorithms  were  coded  using  Symbolic  Manipulation 
methods.  The  work  included  the  analytic  transf ormation  of  the 
"hosted  equations"  (represented  by  the  general  second-order 
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equations  in  3D)  into  general  nonor thogonal  coordinates,  the 
solution  of  these  equations  bv  hopscotch  SQR ,  and  the  coding  and 
solution  of  the  interior  grid  generation  equations  following 
Thompson  et  al .  Symbolic  Man 1  pul  at  1  on  was  used  to  perform  the 
analytical  transformations,  substitute  the  finite  difference 
farms,  gather  and  simplify  terms,  and  write  a  Fortran  subroutine 
for  the  3D  stencil. 

A  rigorous  validation  procedure  was  developed  for  both 
the  hasted  equations  and  the  grid  generation  equations.  This  is 
considered  to  be  a  major  contribution  of  the  research  project. 
The  validation  procedure  consists  in  numerically  testing  the 
truncation  error  convergence.  An  inverse  procedure  was  devised, 
in  which  a  continuum  solution  is  specified,  chosen  so  as  to 
possess  enough  structure  to  exercise  all  the  derivatives  of  the 
operator  and  all  finite  defference  truncation  errors.  The 
non-homogeneous  term  is  chosen  to  give  the  specified  solution. 
By  monitoring  the  truncation  error  as  the  grid  is  successively 
doubled,  all  aspects  of  the  problem  are  verified,  i.e.  the 
tr ansf ormat 1  on ,  the  finite  difference  forms  (validating  the 
second  order  accuracy)  and  the  iterative  solution  procedure.  A 
similar  procedure  is  fallowed  far  the  grid  generation  equations 
by  choosing  a  "solution"  of  a  strongly  twisted  grid.  The  numer¬ 
ical  procedure  is  somewhat  more  complicated  because  the  solution 
for  (x,y,z>  at  each  point  requires  a  nonlinear  transcendental 
solution,  performed  by  Newton-Raphson  iteration. 

These  validation  exercises  also  shed  light  on  a  contro¬ 
versial  question  regarding  the  effect  of  strong  coordinate 
stretching  on  the  order  of  accuracy  of  the  solution  of  the 
hosted  equations.  All  pract i ti oners  had  agreed  that  severe 
inappropriate  stretching  would  increase  the  size  of  the  error, 
but  most  also  were  of  the  opinion  that  the  order  was  also 
degraded,  say  from  second  to  first  order.  Published  analyses 
confirmed  this,  but  the  present  work  contradicted  it,  as  did  our 
straightforward  analysis  in  the  transformed  (logical)  plane. 
The  question  was  settled  under  the  present  contract  by  perform¬ 
ing  systematic  truncation  error  testing  in  ID  using  extreme 
stretching  (involving  even  exponentials  of  exponentials,  requir¬ 
ing  double  precision  calculations).  The  size  of  the  truncation 
error  increased  greatly,  as  expected,  but  it  was  demonstrated 
that  the  order  remained  second  order.  The  other  principals 
involved  have  since  revised  and  clarified  their  analysis,  and 
the  question  has  been  settled. 

The  present  research  effort  also  invalidated  another 
previous  theoretical  (and  intuitive)  result,  that  the  Winslow 
method  (or  homogeneous  Thompson-Thames-Masti n  method)  would 
always  produce  a  nan-folded  grid.  The  present  work  demonstrat¬ 
ed,  by  both  computer  solutions  of  practical  electrode  geometries 
and  by  analysis  of  simple  problems,  that  the  method  can  be  made 
to  fail  in  a  a  predictable  class  of  problems.  The  analysis  has 
bearing  on  the  well  known  difficulties  of  generating  grids  near 
airfoil  trailing  edges,  and  has  been  received  with  considerable 
surprise,  interest  and  incredulity  in  the  grid  generation  commu¬ 
nity. 


A  solution  adaptive  algorithm  and  code  were  developed  tor 
2D  electric  field  calculations.  The  electric  field  calculation 
typical i y  results  in  a  maximum  electric  field  strength  occuring 
on  the  boundary  of  the  region,  e.g.  laser  cavity.  (This  is 
rigorously  true  for  the  vacuum  case,  governed  by  the  simple 
Laplacian  equation. )  Thus,  grid  adaption  does  not  involve 
interior  points  directly  but  rather  the  location  of  grid  points 
on  the  boundary.  This  aspect  of  grid  adaption  is  also  important 
in  fluid  dynamics  but  has  been  largely  neglected.  The  algorithm 
developed  here  uses  a  first  solution  with  boundary  distribution 
determined  by  the  local  surface  curvature.  ( Equi di str i but i on  in 
arc  length  is  the  simplest,  but  a  weighting  in  surface  curvature 
is  also  used.)  With  a  first  estimate  of  the  boundary  electric 
field  values  determined  in  this  grid,  a  second  boundary  distri¬ 
bution  is  chosen  to  adapt  to  this  solution.  The  procedure  can 
be  continued  to  convergence  of  the  adapted  grid,  but  in  practice 
a  single  adapted  grid  is  usually  sufficient.  The  procedure  was 
used  successfully  for  both  steady  state  problems,  and  for 
1 ntra-t 1 me-step  solutions  of  strongly  nonlinear  transient 
pr  ob 1 ems . 

A  family  of  adaptive  rules  was  developed,  with  weighted 
compinations  of  equi di str i but i on ,  boundary  surface  curvature, 
electric  field  strength  (which  is  the  gradient  of  the  solution), 
curvature  of  the  electric  field,  etc.  (Following  recommenda¬ 
tions  by  Eiseman,  the  complete  boundary  surface  curvature  or 
solution  curvature  was  used,  rather  than  just  the  second 
derivative.)  Theoretical  analysis  based  on  Taylor  series  and  a 
criterion  of  minimising  total  error  indicate  that,  with  second 
order  accurate  differencing,  the  grid  points  should  be  distribu¬ 
ted  so  as  to  equi di str i bute  the  (leading  terms  of)  the  trunca¬ 
tion  error.  However,  the  present  work  showed  that  this  is 
impractical,  since  the  resulting  system  is  too  stiff  and  overly 
sensitive  to  noise  in  the  solution.  Several  other  research 
projects  have  come  to  the  same  conclusion,  for  a  variety  of 
different  systems.  Also,  in  the  present  case,  it  was  determined 
that  a  practical  criterion  should  not  even  be  based  entirely  on 
accuracy,  but  should  be  strongly  weighted  to  resolution;  i.e. 
1 arger  error  was  allowed  in  favor  of  more  closely  resolving  the 
maximum  value  of  the  electric  field. 

Although  conceptually  clear,  the  implementation  of  this 
boundary  adaption  involved  a  subtle  procedure  of  "projective 
interpolation" ,  necessitated  by  the  fact  that  the  continuum 
surface  and  the  adapted  boundary  grid  points  do  not  have  the 
same  domain  of  definition.  The  problem  was  successfully  and 
unambiguously  resolved  in  2D,  but  the  analogous  3D  problem  has 
not  been  salved. 

Several  algorithms  were  also  developed  for  the  adaption 
of  the  interior  grid  lines.  The  most  successful  was  a  minor 
extension  of  the  method  of  Thomas  and  Middlecoff.  Unfortunate¬ 
ly,  all  the  algorithms  tested  greatly  slowed  the  iterative 
convergence  of  the  grid  solution  procedure.  This  was  also  true 
of  the  several  algorithms  developed  for  the  interior  boundary 


I 


prob 1 em ,  1 . e . 

i n  ter  1  or  poi nt 
layers  or  inter 
ot  points  is  st 
■:ont  1  nui  t  x  es  l 
angles  unless  t 
boundary.  The 
the  methods  of 
were  only  margi 
was  slow  and  fr 


interior  grid  adaption  with  precise  control  o-f 
s  at  selected  locations  (e.g.  separated  shear 
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nally  successful,  in  that  the  iterative  solution 
aqile. 


Several  iterative  solutions  methods  were  devised  for 
solving  the  grid  generation  equations.  The  GEM  spatial  marching 
methods,  previously  developed  by  the  Principal  Investigator, 
were  successfully  applied  in  semidirect  nonlinear  iterations  to 
the  Thompson  interior  generation  equations  in  2D.  Economies 
result  because  only  one  initial liation  of  the  GEM  algorithm  is 
required  for  the  two  equations  for  x  and  y.  The  method  is  quite 
efficient  for  many  problems,  but  is  fragile,  failing  or  becoming 
very  slow  for  some  geometries  and  for  strong  non-homogeneous 
terms  used  for  interior  grid  control .  Also,  a  previously  unrec¬ 
ognized  deficiency  of  the  marching  methods  was  discovered,  in 
that  gradient  (Neumann)  boundary  conditions  applied  along  bound¬ 
aries  lateral  to  the  march  direction  cannot  be  solved  directly 
unless  the  grid  is  locally  orthogonal  at  these  boundaries. 
Also,  an  improvement  in  the  "smoother"  for  multigrid  algorithms 
for  the  hosted  equations  was  developed  which  is  more  expensive 
of  computer  time  for  simple  problems  but  which  is  not  sensitive 
to  grid  orientation  and  aspect  ratio. 

Initial  progress  has  been  made  on  formulating  an  exten¬ 
sion  of  the  variational  approach  to  grid  generation  pioneered  by 
Brackbill  and  Saltzman.  Symbolic  Manipulation  will  be  utilized 
again  to  handle  the  mathematical  complexity,  especially  in  3D 
and  for  the  non-planar  surface  generation. 

The  Symbol  Manipulation  work  and  the  electric  field 
calculations  were  performed  on  Vax  computers,  and  much  of  the 
grid  generation  work  was  performed  on  a  super-mi crocomputer 
based  on  the  Motorola  68000  chip.  These  combined  systems  have 
proven  to  be  very  effective  for  this  work.  Several  talks  have 
been  given  on  the  experiences  with  these  systems;  interest  has 
generally  been  high. 

The  2D  and  3D  adaptive  grid  generation  algorithms  devel¬ 
oped  in  this  work  have  been  applied  in  the  ELF  codes  (ELectric 
Field)  developed  by  Ecodynamics  under  subcontract  to  Tetra 
Corporation  for  the  Air  Force  Weapons  Laboratory.  These  codes 
are  under  extensive  use  at  AFWL,  at  Tetra  Corporation,  and  at 
university  contractors  for  the  design  and  optimization  of  1 aser 
electrodes  and  switches.  Applications  are  in  the  Pulsed  Power 
area,  including  SDI  research  and  development. 
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University  of  New  Mexico.  Nov.  1984,  Albuquerque,  NM. 

"Computers  vs.  Algorithms",  Patrick  J.  Raache,  Engineering 
Faculty  Seminar,  University  of  Canterbury,  Christchurch,  New 
Zealand,  5  November  1984. 
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'SyntDol  ic  Manipulation  and  Lomputat  l  onal  Fluid  Dynamics11, 
F'atn  cl  J.  Roache,  Engineering  Faculty  Seminar,  University  of 
Canterbury,  Cnr i stchur ch ,  New  Zeal  and ,  7  November  1984. 


(2)  CONSULTATIVE  AND  ADVISORY  FUNCTIONS 

P .  J.  Roache  attended  the  First  Copper  Mountain  Multigrid  Con¬ 
ference  organized  by  Colorado  State  University  on  5-8  April 
1983. 


P.  J.  Koache  and  P.  Eiseman  co-chaired  the  ASME  Mi n i Symposi urn  on 
Grid  Generation  held  at  the  University  of  Houston  on  20-22  June 
1903. 

P.  J.  Roache  and  S.  Steinberg  consulted  with  Dr.  R.  L. 
Rapagnani  and  Capt.  M.  Trout  of  the  Air  Force  Weapons  Laboratory 
on  problems  of  grid  generation  during  1982-1984. 

G.  NEW  DISCOVERIES,  INVENTIONS,  PATENTS,  SPECIFIC  APPLICATIONS 

No  inventions  or  patents  have  resulted  from  this  work. 

The  2D  and  3D  adaptive  grid  generation  algorithms  devel¬ 
oped  in  this  work  have  been  applied  in  the  ELF  codes  (ELectric 
Field)  developed  by  Ecodynamics  under  subcontract  to  Tetra 
Corporation  for  the  Air  Force  Weapons  Laboratory.  These  codes 
are  under  extensive  use  at  AFWL,  at  Tetra  Corporation,  and  at 
university  contractors  for  the  design  and  optimization  of  laser 
electrodes  and  switches.  Applications  are  in  the  Pulsed  Power 
area,  including  SDI  research  and  development. 


The  grid  generation  algorithms  are  also  under 
incorporation  into  a  large  aerodynamics  code  at  AFWL. 
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